How bad metals turn good: spectroscopic signatures of resilient quasiparticles. 
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We investigate transport in strongly correlated metals. Within dynamical mean-field theory, 
we calculate the resistivity, thermopower, optical conductivity and thermodynamic properties of 
a hole-doped Mott insulator. Two well-separated temperature scales are identified: Tfl below 
which Landau Fermi liquid behavior applies, and Tmir above which the resistivity exceeds the 
Mott-Ioffe-Regel value and 'bad-metal' behavior is found. We show that quasiparticle excitations 
remain well-defined above Tpl and dominate transport throughout the intermediate regime Tfl < 
T < Tmir. The lifetime of these 'resilient quasiparticles' is longer for electron-like excitations, and 
this pronounced particle-hole asymmetry has important consequences for the thermopower. The 
crossover into the bad-metal regime corresponds to the disappearance of these excitations, and has 
clear signatures in optical spectroscopy. 



The transport properties of metals with strong electron 
correlations are unconventional and poorly understood 
theoretically. Two facts regarding the temperature de- 
pendence of the resistivity are frequently observed, (i) 
Fermi-liquid (FL) behavior p oc T 2 only holds below a 
temperature Tfl which is low compared to bare elec- 
tronic energy scales, (ii) At high temperatures the re- 
sistivity is large and reaches values exceeding the Mott- 
Ioffe-Regel (MIR) limit. This 'bad-metallic' behavior [1] 
signals the breakdown of a quasiparticle (QP) description 
of transport, since the associated mean- free path I would 
be smaller than the lattice spacing. This is observed in 
many materials: Sr2Ru04 has Tfl — 20K, while the MIR 
value is reached at Tmir — 800K (using fcpZ ~ 1 as the 
MIR criterion); in LiV 2 4 , T F l is a few degrees Kelvin, 
while Tmir is several hundreds, etc. (see [H [3] for re- 
views). In two-dimensional organic materials, Tfl and 
Tmir are closer but still distinct scales [3HS]- 

These observations raise the following questions. Why 
is Tfl much lower than Tmir and what determines its 
value? Up to which temperature do QPs exist and what 
are the signatures of their disappearance? And, most 
importantly: how should one think of transport not only 
in the bad-metal, but also in the intermediate regime 
?fl ;$ T < Tmir where the resistivity does not fol- 
low Landau's T 2 behavior, but is still smaller than the 
MIR value? These questions also apply to cuprate su- 
perconductors, where the observation of quantum oscil- 
lations [7J [5] and T 2 behavior in transport [jH HH] and 
optics [TT] have rejuvenated the relevance of FL states 
with a low Tfl (possibly with angular dependence along 
the Fermi surface). 

In this article, we answer these questions in a partic- 
ularly simple setting: a hole-doped Mott insulator de- 
scribed with dynamical mean- field theory (DMFT). Our 



most striking finding is that well-defined QP excitations 
survive well above the range of validity of FL theory. 
For over a decade above Tfl the transport can be ac- 
curately described in terms of these 'resilient quasiparti- 
cles' (RQPs). The high-temperature MIR crossover into 
the bad-metal regime is associated with their gradual ex- 
tinction, which has a clear signature both in the single- 
particle spectral function and in optical spectroscopy. 
In a hole-doped Mott insulator, the RQPs come with a 
strong particle-hole asymmetry: electron-like excitations 
are longer-lived than hole-like ones. This has direct con- 
sequences for the thermopower. 

Previous DMFT work has investigated transport [T21 - 
ITB] and optical conductivity [T7J, but not the precise 
temperature-dependence of the self-energy and of the 
momentum-resolved spectral function which reveals this 
intermediate RQP regime. We note that in the half- 
filled case relevant to organic compounds 4, 5J, the high- 
temperature state is insulating-like, and hence the tem- 
perature window where this regime can be seen is nar- 
rower. 

We solved the DMFT equations [18 for the hole- 
doped Hubbard model using highly accurate numerical- 
renormalization group (NRG) [TU [2D] and continuous- 
time quantum Monte Carlo [21] [22] techniques. In 
DMFT the real part of the optical conductivity reads 

<t(w) = 2n ^f du/F UtU > J de<Z>(e)A k (uj')A k (u>'+Lj) (1) 

where F UtU i — [/(&/) — f(u) + u>')]/w with f(oj) the Fermi 
function. A^(lu) = — (l/7r)Im[Lj + /i — £k~ S(cj)] -1 is the 
single-particle spectral function, with the energy of the 
state in the band, the retarded self energy and fi the 
chemical potential. <J>(e) = (1/V) J2k(d £ k/dk x ) 2 5(e - 
£k) contains the information about velocities. We used 
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FIG. 1. (a) Temperature-dependence of the resistivity for 
several doping levels 8 (U/D — 4, as in all figures). The 
MIR value as denned in the text is reached at a temperature 
Tmir indicated by plain arrows. Inset: resistivity at low tem- 
peratures vs. T/8D revealing the T 2 behavior (dashed line) 
below Tfl (empty arrow), (b) Determination of Tfl by a 
scaling plot of Z(T)/Z{T -> 0) vs. T/SD. Here, Z(T)" 1 = 
1 - dReE(uj,T)/doj. (c) The different regimes: FL (blue) 
for T < Tfl, bad metal (red) and intermediate RQP regime. 
The crossover into the bad metal is gradual: the onset of red 
shading corresponds to the optical spectroscopy signatures 
discussed in the text, while the red points indicate where pmir 
is reached. The thin dashed line indicates the 'knee' in p(T). 



a semicircular density-of-states (DOS) with a half-width 
D, and the corresponding sum-rule preserving expression 
$(e) = $(0) [l - (e/D) 2 ] 3/2 . In the following, resistivity 
will be expressed in units of the MIR value defined as 
1/pmir = e 2 &(0)/HD. This is consistent with the crite- 
rion kp I = 1 for a parabolic band in two dimensions, for 
which the conductivity a — (kpl) e 2 $(ef)/feF. 

Fig. [l£i displays the resistivity vs. T for several dop- 
ing levels, and Fig. [T]; summarizes the three different 
regimes. The high temperature crossover into the 'bad 
metal' regime is a gradual one. Its onset can be defined 
from changes in the optical conductivity (see below), 
while the temperature Tmir at which the MIR value is 
reached (as defined above) provides an upper estimate. 
The resistivity increases linearly (with a positive inter- 
cept) in this high-T regime, with no evidence of satura- 



tion. Our data reveal that the onset of the bad metal 
regime occurs at the Brinkman-Rice scale (~ SD), which 
is a measure of the kinetic energy of QPs (and hence of 
the quantum degeneracy scale). Note that, within local 
DMFT, both the QP weight and mass renormalization 
are proportional to doping Z — m/m* ~ S. 

At very low T < T F l, FL behavior is found. We 
determined Tfl from several complementary criteria: 
the onset of T 2 resistivity (inset of Fig. [l^,), scaling of 
Im£(u;, T)/T 2 vs. w/T (see on-line supplementary |23|). 
scaling of Z(T = 0)/Z(T) vs. T/SD (Fig.[l]b). We find 
Tfl — 0.05SD, also proportional to doping, but with a 
much smaller prefactor than Tmir - Correspondingly, the 
resistivity at Tfl is way smaller than the MIR value. A 
low-T expansion yields: p(T)/p M m ~ 6.3(T/<5T>) 2 + • • • 
and hence /c(Tfl)/pmir — 0.016. 

There is thus a wide temperature range XpL T < 
Tmir in which transport does not follow the T 2 FL behav- 
ior, although resistivity is still substantially smaller than 
Pmir- Right above Tfl, the resistivity increases approx- 
imately linearly (with a negative intercept). A knee- like 
feature is observed at a temperature T (dashed line in 
Fig. [TJ:) , above which the high-T regime gradually sets in, 
in which p(T) has linear temperature dependence (with 
a positive intercept), as can be shown from a high-T ex- 
pansion p~5| ri6] . 

The above findings raise the following question: what 
are the charge carriers in the intermediate metallic regime 
Tfl < T < Tmir? We depict in Fig. [2] the momentum- 
resolved spectral function A-^(lu) at selected tempera- 
tures as energy distribution curves {A^{uj) vs. u> for differ- 
ent £k's). Also shown is the momentum-integrated DOS. 
These results reveal a remarkable fact: well-defined QP 
excitations exist throughout this intermediate regime, 
way above the FL scale. Our definition of the term 'quasi- 
particle' is a pragmatic one: we mean that A^uj) displays 
a well-resolved peak in the vicinity of the Fermi level, in 
addition to a lower Hubbard band (LHB) and an upper 
Hubbard band (UHB). 

For T < T FL (e.g. T/D = 0.0025 in Fig. [2} sharp peaks 
are seen close to the Fermi energy (w = 0), characteristic 
of long-lived Landau QPs. ForT > T FL (e.g. T/D = 0.05 
and T/D = 0.2 in Fig. [2]), the peaks broaden and the 
RQPs are visible mostly on the unoccupied -side of the 
FS [u > 0). Finally, for T > T MIR , all QP features are 
gone, and only Hubbard satellites remain in the spectra. 

Examination of the self-energy (Fig. [3]) helps under- 
standing the nature of the QP excitations, as well as of 
the different transport regimes. In local DMFT, — ZImS 
can be interpreted at low-w as the inverse of the QP life- 
time, and — ImE as the transport scattering rate. Fig.[3]b 
displays Im£(w,T) vs.T for different excitation energies 
uj. At uj = (thick curve), FL behavior ImE ~ T 2 applies 
at low-T, corresponding to very long-lived QPs. Note 
that strict FL behavior breaks down already below the 
temperature at which -ZImE(0,T) ~ T [21]. At finite 
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FIG. 2. (a) Temperature evolution of the total DOS for 
S — 0.2. (b) Momentum- resolved spectral functions. The 
shaded area [SksT fiksT] indicates the states with a signifi- 
cant contribution to transport. 




FIG. 3. Self-energy and particle-hole asymmetry, (a) ImE(u) 
for different temperatures, (b) Temperature dependence of 
ImS(wc) for io c = —0.5, —0.4, . . . , —0.1 (turquoise), u) c = 0.0 
(thick black) and uj c = 0.1,0.2,...,0.5 (green). Below the 
dashed line Z(T)Im£ < T. (c) ReE(cj) at T/D = 0.0025, 
and the two 'kinks' (arrows). 



frequency, the hole-like excitations have higher scattering 
rate than electron-like ones (Fig. [3]a,b). At ~ 0.08Z?, 
the curves of ImS vs.T for different positive w's display 
a crossing point. Above T*, the scattering rate is a de- 
creasing function of frequency: low-energy electron-like 
excitations with finite uj > have a longer lifetime than 
those at u = 0. These finite-energy uj > excitations 
provide the largest contribution to conductivity in the 
intermediate RQP regime [35] • Their inverse lifetime de- 
pends weakly on temperature for T > T* (plateau-like 
behavior in Fig. [3]b) and remains much smaller than the 
bandwidth and at most comparable to fc^T. For early 
considerations on a QP description of transport beyond 
the FL regime in the context of electron-phonon interac- 
tions, see [25]. 

The strong electron-hole asymmetry has also other in- 
teresting consequences. Because in the RQP regime the 
uj < states are strongly damped, the Fermi surface 
as determined by the maximum intensity of A^{ui = 0) 
'inflates' to a larger volume than the T — Luttinger 
volume [53]. From Fig. [3]c, it is also seen that the 
deviation from linearity of ReE (~ So + w(l — l/Z) 
at low uj) defines two distinct energy scales, — w_ for 
hole-like and uj + for electron-like excitations, leading to 
'kinks' in QP dispersions, as documented by previous 
studies [271 125] . We note that the smallest kink energy 
lu + <C w_ sets the scale for deviations from FL behavior 
(u>-\- ~ 7tTfl)- Using quite different theoretical methods, 
previous studies [29l [30] have also emphasized the im- 
portance of particle-hole asymmetry in hole-doped Mott 
insulators. 

A sensitive probe of the particle-hole asymmetry is the 
Seebeck coefficient (thermopower) Q(T) shown on Fig. 51. 
Strikingly, the sublcading particle-hole asymmetric terms 
in the low-frequency expansion of ImS modify the slope 
of Q(T) at low-T by a factor of about two, as compared 
to a naive FL theory estimate (thick dashed line) that 



would only retain terms 



+ (ttT) 2 . This effect was 



anticipated in Ref. 31J and is shown here to be quantita- 
tively important. The plateau behavior of the scattering 
rate of the RQPs, discussed above, is also responsible 
for Q(T) still increasing in an electron-like manner up 
to T ~ T*. At a higher temperature within the RQP 
regime Q(T) changes sign and, when entering the bad- 
metal regime, approaches the simple Heikes estimate for 
D<T < [/ (Fig. The atomic Kelvin formula [H[32] 
successfully describes the thermopower there, which can 
thus be taken as another fingerprint of a bad-metal. The 
accuracy of approximate formulas for thermopower has 
been tested also in other studies [33] [31] . 

It is interesting to observe how the different trans- 
port regimes relate to thermodynamic observables. On 
Fig. |4]wc display the entropy S(T), the kinetic energy 
K(T), and the Curie constant Txi oc associated with the 
local (q-integrated) magnetic susceptibility. The entropy 
as well as the Curie constant reach remarkably high val- 
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FIG. 4. Entropy (in units of Jcb), Txioc (Ms = l)i aQ d See- 
beck coefficient (in units of fes/|e|), for S = 20%. Horizontal 
dashed lines indicate the atomic (Heikes) limit for the entropy 
and thermopower. Vertical dashed lines denote Tfl and Tmir. 
The thermopower at low-T is compared to the FL estimate 
(thick dashed) in which the particle-hole asymmetry of ImE 
is neglected. Kinetic energy in units of D (right scale). Inset: 
Entropy and kinetic energy versus temperature on a log scale. 




FIG. 5. Optical conductivity for 5 — 20%. (Inset:) Optical 
spectral weight integrated up to fl, normalized to the kinetic 
energy. 



ues already below Tmir- In the RQP regime, the sys- 
tem thus has to be thought of in terms of two fluids: a 
mixture of local moments and of the resilient QP states. 
Above Tmir, the entropic contribution to the free energy 
overcomes the kinetic energy gain. The system is fully 
incoherent, its entropy approaches the atomic limit and 
the Curie constant saturates. 

Optical spectroscopy (Fig. [5]) can be used to detect the 
two crossovers, between the FL and the RQP regime and 
from the latter into the bad-metal. For T < Tfl, <t ( uj ) 
displays a narrow low-frequency peak which decays as 
1 juJ 1 . This Drude peak corresponds to optical transitions 



involving only QP states and has a spectral weight pro- 
portional to doping level S. In the FL regime, the Drude 
peak is well separated from a higher frequency 'hump' 
(at ui/D ~ 0.5 in Fig. [5]) which corresponds to transi- 
tions between the LHB and QP states [T7] . The distance 
of the LHB to the Fermi level, of order /i (a fraction 
of the bandwidth) sets the energy scale for these tran- 
sitions. This typically corresponds to the mid-infrared 
range in narrow-band correlated materials. The crossover 
out of the FL regime into the 'resilient metal' regime 
leads to a broader low-frequency peak whose frequency 
dependence is no longer 1/u; 2 . (For T/D = 0.02 and 
0.02 < lo/D < 0.2, it can be fit to l/oj a with a w 1.2; 
We also notice that at low-T, <t(uj) has a non-Drude 'foot' 
in this frequency range.) Above T* the Drude and mid- 
infrared features merge. 

For a rather extended temperature range into the RQP 
regime, the spectral weight redistribution upon increas- 
ing T takes place essentially entirely between the low- 
energy QP states and the 'mid-infrared' feature, see inset 
of Fig. [5]. For T > T* (warmer part of the intermediate 
regime) , some spectral weight transfer to higher frequen- 
cies starts taking place as well. The crossover into the 
bad-metal regime at Tmir is signaled by two changes in 
a(w). First, the isosbestic crossing point (at 00/ D ~ 0.5) 
is lost and the low-frequency non-Drude peak is replaced 
by a very broad peak. This peak corresponds to tran- 
sitions involving only the LHB. Furthermore, spectral 
weight is now redistributed over a considerable energy 
range, extending all the way to the UHB (inset). In the 
bad-metal, the kinetic energy (to which <j(oj) integrates) 
is strongly dependent on temperature (Fig. [4]) and the 
broad peak has correspondingly a height which contin- 
ues diminishing with T, which implies that the resistivity 
does not saturate [2]. 

In summary, our study reveals that resilient QP ex- 
citations persist well above Tfl- They control trans- 
port properties until they disappear at T ~ Tmir, a s 
revealed e.g. by optical spectroscopy. The coexistence of 
QP states with localized magnetic moments which carry 
large entropy is intriguing, and demands closer theoreti- 
cal investigation including antiferromagnetic correlations 
beyond DMFT. For hole-doped Mott insulators, a pro- 
nounced particle-hole asymmetry is found. This calls 
for new momentum-resolved spectroscopic probes which 
would be able to access the 'dark side' of the Fermi sur- 
face and for closer investigations of the electron-doped 
materials where the signatures of this regime could be 
seen using conventional ARPES. 
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Supplementary information 



Self-energy scaling 

In a Fermi liquid, Im£(w,T) oc [u 2 + (irT) 2 ] and 
therefore the self-energy obeys the following scaling in 
uj/T: 



1000 



L»ImS 



A [it 2 + (uj/T) 2 ] . 



(A.2) 



In Fig. [fjj we show — Iml](w, T) /T 2 as a function of cj/T 
for different temperatures. As the temperature is lowered 
below Tfl ~ 0.01D, the curves collapse on a parabola, 
confirming the expected Fermi-liquid scaling law. Find- 
ing such a scaling is actually a stringent test on the nu- 
merical data. Very precise quantum Monte Carlo data 
(analytically continued with Pade approximants) and a 
strict control of the chemical potential were needed to 
obtain these results. 

As the temperature is raised, the positive-frequency 
side quickly deviates from the scaling form, revealing a 
particle-hole asymmetry already for temperatures above 
T/D = 0.005. The deviations appear at a scale ui + 
7tTfl. On the negative- frequency side the self-energy fol- 
lows the quadratic behavior much more robustly and de- 
viates from the scaling function only at about T/D = 
0.02. Notice that the corrections on the uj > side grow 
linearly with temperature. This leads to a Seebeck coef- 
ficient which is linear in T at low temperatures, but with 
an enhanced slope as compared to the result one would 
get if these corrections were neglected. 

The transport probes an energy window of a few (say 
from -5 to 5) ksT . In this energy window the self-energy 
starts to deviate appreciably at T/D — 0.01. This is 
where the resistivity (within the precision of our data) 
visually departs from the T 2 law. 



Temperature evolution of momentum-resolved 
spectra 




FIG. 6. The self energy scaling as uj/T for the doped Hubbard 
model considered in the main text (Bethe lattice, U = AD, 
dop ing 5 = 0.2). At low temperatures self energy follows 



Eq. 



A.2 



A is found to be close to Z 



(For doping 



0.2, Z = 0.22.) 




FIG. 7. Contour map of momentum-resolved spectra A(ek, ui) 
at various temperatures for the same parameters as in Fig[6] 



In Fig. [7] we plot the temperature evolution of the 
momentum-resolved spectra using a color-map where 
bright (dark) colors indicate high (low) values of Aj.{u). 
At a low temperature, the data displays three peaks cor- 
responding to: the lower Hubbard band (LHB) at energy 
wlhb ~ — /•*: a quasiparticle peak (QP) in the vicinity 
of uj ~ 0, and upper Hubbard band (UHB) at energy 
wtjhb ~ U — fi. On the top left panel one can clearly 
distinguish the quasiparticle band from the lower Hub- 
bard band (the upper Hubbard band is above the upper 
energy bound of the plot). At very small frequencies, the 
slope of the dispersion is found to be approximately 5 
times smaller than the bare one, as dictated by2w 0.2. 
On the negative-frequency side this holds almost until 



the bottom of the band. On the positive-frequency side, 
instead, the kink at ui + is rapidly encountered. Above 
this kink, the slope of the band dispersion increases to 
about half the bare dispersion slope. 

As the temperature is increased, the quasiparticle band 
broadens but remains clearly visible (left panels). On 
approaching Tmir the quasiparticle part of the spectra 
progressively merges (right top two panels) into the lower 
Hubbard band until (right bottom two panels) no signal 
on top of the Hubbard band can be discerned. 

Note that the quasiparticle band crosses the Fermi en- 
ergy at different momenta as the temperature increases. 
Identifying the Fermi surface with the momenta at which 
the spectral intensity of the quasiparticle band is maxi- 
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FIG. 8. The evolution of momentum-distribution curves at 



the Fermi level Eq. (A.3l 



mal leads to a conclusion that the Fermi volume inflates 
as the temperature is increased. Note that the number 
of particles is fixed so that with this identification of the 
Fermi surface the Luttinger theorem is only obeyed at 
very low temperatures. 

To elaborate on this, we plot on Fig.[8]the momentum- 
distribution curve at u) = 0, 



signature of well-distinguishable resilient quasiparticles. 



Thermopower at high temperatures and comparison 
to approximate formulas 

In Fig. [lO] we show the Seebeck coefficient over a larger 
temperature window. The Seebeck coefficient calculated 




0.4 0.6 
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FIG. 9. Temperature dependence of fj, e s = n(T) — Re£(aj = 
0,T). The chemical potential for non-interacting chemical 
potential is shown by a dashed line. 



A(£fe ' 0) 7T (n - Rc£(0, T) - e k y + (lm£(0, T)) 2 

(A.3) 

for several temperatures. This spectral function has the 
shape of a Lorentzian centered at /i — Re£(0, T) — with 
a width ImS(0,T). At very low temperatures a sharp 
peak lies at the non-interacting chemical potential, which 
fulfills the Luttinger theorem. When the temperature 
increases, the peak moves to higher momenta. 

An alternative way to track this change is to look at 
the renormalized chemical potential (Lt e g = fi— Re£(0,T) 
as a function of the temperature as shown in Fig. [9] 
In the Fermi-liquid regime, n e g essentially follows the 
non-interacting chemical potential (shown with a dashed 
line). At elevatated temperatures /j, e s rapidly increases. 

An important lesson here is that observing a well- 
defined quasiparticle peak does not imply that the system 
has reached the Fermi liquid regime. As Fig. [9] shows, 
the chemical potential in this intermediate-temperature 
metal can be quite far from the Fermi energy, despite a 



using the Kubo formula is plotted with a thick line 
and compared to various estimates. The Kelvin formula 
dfi/dT (using /i calculated within DMFT) overestimates 
the magnitude of the thermopower in the low-T regime 
but is a good approximation of the exact result above 
T* w 0.08D. For comparison, we also plot the corre- 
sponding atomic estimate using the Kelvin formula, but 
using the /z obtained in the atomic limit. We note that 
the two expressions essentially match above T/D = 1. 
Note that the Kubo result starts to deviate significantly 
from the atomic estimate only when entering the resilient 
quasiparticle regime. 

Finally, we plot the Heikes estimates. The results ob- 
tained from the DMFT chemical potential using NRG 
and continuous-time interaction expansion Monte Carlo 
(CTINT) are plotted (full green line and symbols). The 
Heikes formula is found to approximate the thermopower 
worse than the Kelvin formula. For comparison, also the 
atomic Heikes estimates (thick dashed) as well as the 
asymptotic U — > and U — > oo Heikes values (horizontal 
lines) are shown. 
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FIG. 10. Seebeck coefficient calculated using the exact 
Kubo formula (thick black line) compared to the approxi- 
mate Kelvin formula (red line) and the Heikes formula (green 
line and green symbols). The atomic Kelvin estimate (dot- 
ted) and Heikes estimate (thick dashed) interpolate between 
the asymptotic U — > oo and U — > Heikes values. 



